Astronomy &: Astrophysics manuscript no. parizotEtAlXrayRims February 5, 2008 

(DOI: will be inserted by hand later) 



Observational constraints on energetic particle diffusion in young 
SNRs: amplified magnetic field and maximum energy 

E. Parizot\ A. Marcowith^, J. Ballet^, and Y. A. Gallant*^ 



\^ [ ^ IPN Orsay, CNRS/Universite Paris-Sud, 91406 Orsay Cedex, France 

■ ^ Centre d'Etude Spatiale des Rayonnements, 9 av. du Colonel Roche, 31028 Toulouse Cedex 4, France 

O ' ^ DSM/DAPNIA/SAp, CEA Saclay 91191 Gif sur Yvette Cedex, France 

CN ! ^ LPTA, CNRS/Universite Montpellier II, 34095 Montpellier Cedex 5, France 



(N 



Received ??; Accepted ?? 

Abstract. Constraints on the diffusion and acceleration parameters in five young supernova remnants (SNRs) 
are derived from the observed thickness of their X-ray rims, as limited by the synchrotron losses of the highest 
energy electrons, assuming uniform and isotropic turbulence. From a joint study of the electrons diffusion and 
advection in the downstream medium of the shock, it is shown that the magnetic field must be amplified up to 
^ . values between 250 and 500 fiG in the case of Cas A, Kepler, and Tycho, or ~ 100 fiG in the case of SN 1006 

^Nf^ ' and G347.3-0.5. The diffusion coefficient at the highest electron energy can also be derived from the data, by 

^SJ , relating the X-ray energy cutoff to the acceleration timescale. Values typically between 1 and 10 times the Bohm 

■ diffusion coefficient are found to be required. We also find interesting constraints on the energy dependence of 

CO ' the diffusion coefficient, by requiring that the diffusion coefficient at the maximum proton energy be not smaller 

, than the Bohm value in the amplified field. This favours diffusion regime between the Kraichnan and the Bohm 

regime, and rejects turbulence spectrum indices larger than ~ 3/2. Finally, the maximum energy of the accelerated 
particles is found to lay between 10^^ and 5 10^'' eV for electrons, and around Z x 8 lO^** eV at most for nuclei (or 
. ~ 2.5 times less if a Bohm diffusion regime is assumed), roughly independently of the compression ratio assumed 

CL(' at the shock. Even by taking advantage of the uncertainties on the measured parameters, it appears very difficult 

for the considered SNRs in their current stage of evolution to produce protons up to the knee of the cosmic-ray 
spectrum, at ~ 3 10^^ eV, and essentially impossible to accelerate Fe nuclei up to either the ankle at ~ 3 x 10^* eV 
or the second knee at 



o 



Key words. Acceleration of particles - Magnetic fields - Cosmic rays - ISM: supernova remnants 



H 1- Introduction the resonant interaction with MHD waves. The value and 

- - structure of the inhomogeneous magnetic field is thus cru- 

High-angular resolution X-ray observations of young su- cial to determine D{E). In the absence of a fully predic- 

pernova remnants (SNRs) show very thin rims of emis- tive theory of wave generation and particle diffusion inside 

sion, associated with the forward shock of the supernova SNRs, it is generally assumed that the diffusion coefficient 

(SN) expanding in the interstellar medium (ISM). This is close to the so-called Bohm limit at all energies, defined 

emission is most probably rel ated to the syn chrotron emis- by D^{E) — rLw/3, where tl = p/qB is the Larmor ra- 

sion of high-energy electrons (iBallet 1 1200,4 and references dius of the particle of momentum p in the field B, and q 

therein) accelerated at the shock by the well-known diffu- and v are its charge and velocity. This value of D corre- 

sive shock acceleration (DSA) mechanism, which requires sponds to a mean free path of the charged particles equal 

the scattering of electrons back and forth through the to the Larmor radius, which is thought to be the lowest 

velocity discontinuity associated with the shock. While possible value for isotropic turbulence. However, particles 

this mechanism is well understood qualitatively and has of different energies resonate with field fiuctuations at dif- 

been extensively confronted with multi-wavelength obser- ferent scales, and the magnetic field intensity to be used 

vational data, as well as direct measurement in interplan- to calculate tl (and thus Db) should be an effective value 

etary shocks, very little is known about one of its most at the resonant scale (A ~ tl). The diffusion coefficient 

fundamental ingredients: the particle diffusion coefficient, therefore depends in principle on the power spectrum of 

The diffusion mechanism is related to the deflection ^^^e magnetic field. In particular, D{E) oc E^'^ is expected 

of charged particles in the ambient magnetic field and ^^e quasi-linear regime in the case of a Kolmogorov- 
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like spectrum of turbulence. Recent numerical analyses 
of particle diffusion in isotropic turbulent magnetic fields 
have confirmed analytical results in their range of valid- 
ity a nd shown that the Bohm regime does not generally 
hold l|Casse et al. Il2002t IParizot ll20oi ICandia fc Roulet I 
120041) . even in the strong turbulence limit, and can be ap- 
proached only at the critical energy where tl = Ac, the 
coherence length of the field. The question thus remains 
whether the Bohm approximation is relevant to the case 
of DSA in SNRs. 

A striking property of the DSA mechanism, however, 
is that in the steady-state regime and the so-called test- 
particle approximation, the global shape of the energetic 
particle (EP) energy spectrum does not depend on the 
exact value of the diffusion coefficient, nor even on its 
variation with both energy and space. Many studies could 
thus be made without precise knowledge of the diffusion 
mechanism. In particular, the 'universal' spectrum 
could be obtained in the test-particle limit, and some in- 
sights were gained into the non-linear regime where the 
EPs influence the shock structure (see e.g. the review by 
Jones and Ellison, 1991). The acceleration rate, however, 
is directly related to the spatial diffusion coefficient: the 
smaller this parameter, D{E), the smaller the acceleration 
time, Tacc- To some factor of order unity, race — D{E)/V2i, 
where V^h is the shock velocity. For this reason, the value 
(and energy dependence) of D{E) governs the highest en- 
ergy which can be reached in a SNR. 

Another aspect of diffusion, of course, is the trans- 
port of particles in space! By allowing one to draw 
maps of the non thermal emission, high-resolution ob- 
servations of SNRs can provide information about the 
EP t ransport properties ijvan der Swaluw fc Achterberg I 
I2OO4 . The other ingredients governing particle transport 
are advection and energy losses. In the case of electrons, 
the latter are also related to the ambient magnetic field, 
and thus are not independent of the diffusion coefficient. 



In this paper, we extend pr evious work on spe- 
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Yamazaki ct al. 
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constraints that high-angular resolution observations can 
set on two key parameters: the magnetic field and the dif- 
fusion coefficient. We discuss the interplay between the 
various phenomena governing electron transport in SNRs, 
assuming that the turbulence is isotropic: investigation of 
anisotropic effects are beyond the scope of the paper. We 
make the further assumption that the X-ray rims thickness 
is limited by the radiative losses and leave for future work 
the discussion of the impact of the turbulence relaxation 
downstream of the forward shock on the X-ray filaments 
morphology (see Pohl et al. 2005). We apply our analy- 
ses to five young SNRs observed by Chandra and discuss 
the resulting maximum energy of electrons and protons 
in SNRs. For all numerical applications, we assume the 
reference values summarized in Tabled 



2. Electron transport at the SNR shock 

2.1. The size of the X-ray rims 

In the DSA theory, energetic electrons are assumed to be 
scattered by magnetic inhomogeneities (or MHD waves) 
hosted by the underlying plasma, and diffuse roughly 
isotropically in the plasma rest frame. With respect to 
the shock front, two processes compete: i) global advection 
with the fluid, towards the shock upstream and away from 
it downstream, and ii) spatial diffusion which allows a frac- 
tion of the particles advected downstream to nevertheless 
cross back the shock and gain energy through the stan- 
dard first order Fermi mechanism. Diffusion also makes 
the upstream particles explore a region of size D/Vgh 
ahead of the shock, where Vsh is the shock velocity. In ad- 
dition, the highest energy electrons suffer from efficient 
energy losses due to synchrotron emission in the local 
magnetic field, which limit their maximum energy and 
lead to an assumed exponential cut-off in their spectrum: 
fc{E) cx £'"'•' X exp(— £^/£^c,max)- Because of these losses, 
the highest energy electrons cannot travel far from the 
acceleration region and their emission must be limited to 
very thin regions just behind the shock, seen as sharp and 
narrow rims in X-ray. 

The observed thickness of the rims, Ai?obs, is related 
to their actual size by a projection factor, P, which de- 
pends on the geometry of emission region. From simple 
arguments based on the transport equation of the accel- 
erated electrons (see below), one can infer that the elec- 
tron distribution downstream has an exponential profile: 
fcix) oc exp(— x/Ai?i-i,ii), where x is the radial coordi- 
nate measured from the shock front, with increasing values 
downstream, and Ai?rim defines the characteristic scale of 
the electron distribution. Note that Ai?i.im depends on 
the electron energy. Here, we always refer to the observed 
rims at E^ = 5 keV, corresponding to electrons of energy 

E = £'o,obs. 

In the ideal case of a spherical shock, an exponen- 
tial emission profile translates into a projected thickness 
Ai?obs — 4.6Ai?rim, where Ai?obs is defined as the size 
of the region where the observed brightness is larger than 
half of the maximum (Ballet, 2005). This was obtained 
by numerical calculation in the limit Ai?obs ^ -Rsh, the 
SNR shock radius. Note that we neglected the contribu- 
tion of the upstream emission, which should indeed be 
much thinner or even inexistent at the X-ray energies con- 
sidered, because of the lower magnetic field upstream (see 
below). In the following, we adopt the above value of the 
projection factor, P = 4.6, and use the reduced projection 
factor P = P/4.6 to track how the results depend on the 
assumed geometry. 

2.2. Synchrotron loss time 

Electrons of energy E = jrricC^ propagating in a magnetic 
field B with a pitch angle a emit a synchrotron radiation 
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Name 


Distance 


Shock speed 


Projected width 


Age 


Cut-off energy 


Cas A 
Kepler 
Tycho 
SN 1006 
G347.3-0.5 


3.4 kpc ' 
4.8 kpc ^ 
2.3 kpc ^ 

2.2 kpc 12 

1.3 kpc i« 


5200 km/s ^ 
5400 km/s ^ 
4600 km/s ^° 
2900 km/s 
4000 km/s ? 


0.05 pc (3") ^ 
0.07 pc (3") ^ 
0.05 pc (4") ^ 
0.2 pc (20") " 
0.25 pc (40") 


320 yr ? * 
400 yr 
430 yr 
1000 yr 
1620 yr ? 


1200 eV ^ (whole SNR) 
900 eV * 
290 eV " 
3000 eV 1^ 
2600 eV 1^ 



Table 1. Characteristics of the non-thermal emission behind the blast wa ve in young SNRs (s ee Ballet, 2005, fo r 
f urther detail). The numerical values are taken (or directly inferred) from: ^ iReed et al l lll995l) . ^ IVink et all lIlQQal) 
^ IVink fc Laming I ll20 031 from a fit to the spectrum of the whole SNR, ^ lAshworthl lIlQSOl). ^ iRevnoso fc Goss I 
l)l999(l . '' iHughes I lll99 91)+ distance, ^ De cour chelle (2003) (p rivat e communicat i on) , ^ ICassam-Chenaf et al. I (l2004al) 
in the southeast ^ ISmith et a,1 1 ()l99l|^ . ^° iHughesI ll20f)(|. " iHwa.ng et aJ I (|2n02l). IWinkler et al I tOQc$\. 
IChavamian et al I ll2002l) in the northwest,^ ^'^ Bamba et al I ((200,'^ in the northeast, iR.othenflug et al I ll2 004) in the 
northeast. ""^^ ICassam-Chenai et al" * ( 2004b) , from v — Xr/t with the expansion parameter A = 2/3. ^^ iLazendic et al I 



an^ei 

— t ; , 

in the northwest, IWang et al . (,,1997, 



at a typical energy of 



3/i 2 6^ 
4-;— 57 — sma, 
47rp rric 



(1) 



where ^ ~ 0.29 (e.g. Longair, 1994). This can be inverted 
to give the typical energy of an electron emitting photons 
at energy = i?-y,keV keV. We shall assume here an ap- 
proximate "mono-energetic" (i.e. one-to-one) conversion, 
as given by: 



E, 



(2) 



where Bioo is the magnetic field in units of 100 fiG (The 
impact of this simplification on our numerical results is 
negligible at the level of precision of the present approach) . 
For the observed X-ray energy, E^ = 5 keV, we thus have: 



Se,obs=^(49TeV)xSio;/^ 



(3) 

When the electron leaves the acceleration process (i.e. 
it will not cross the shock front again), it cannot radi- 
ate at energy E-y for a time longer than the synchrotron 
energy loss time, Tgyn, due to the synchrotron radiation 
itself (note that synchrotron losses are usually expected 
to be larger than inverse Compton losses, because of the 
large magnetic field). The corresponding energy loss rate 



is written {dE/dt)s 



\aTC X emagn X , whcrc 



Cmagn = -B~/2/io is the magnetic energy density, ctt is the 
Thomson cross section and 13 = v / c ^ 1 in the case of 
interest here. Thus, noting -EtcV the electron energy in 
units of 10^2 gV: 



E 



(1.25xl03yr)x£;TeVx5roo, (4) 



- (d£;/dt)sy„ 
or, at the energy corresponding to the observed rims: 



'syn 



(Se.obs) ^ (25 yr) x B^^^^ . 



(5) 



2.3. Advection downstream 



where r is the shock compression ratio (r = 4 for an un- 
modified strong shock, but it can be larger if non-linear 
effects are important, e.g. Ellison et al., 2004). The size of 
the advection region, where the electrons emit at 5 keV, 
is thus simply given by Ai?adv = Ki x Tsyn(i?c,obs)- It can 
also be shown more formally that this is indeed the scale 
of the rims downstream, Ai?i.i,„, obtained when neglecting 
electron diffusion. The standard transport equation at a 
parallel shock is written (Volk ct al., 1981): 



f 



dt 



dx 3 dp dx dx dx ' 



+ g, (6) 



' syn 



where the last term, g, is an unspecified source term, and 
the loss term on the right hand side has been given a sim- 
plified form, corresponding to random catastrophic losses 
after a lifetime Tgyn. The stationary version of Eq. © 
in the downstream region and without diffusion is sim- 
ply Vdidf/dx) = — //Tsyn, with the obvious solution 
fe{x) (X exp(-a:/ydrsy„). 

Taking Tgyn from Eq. jSJ, one obtains: 



Ai?adv ^ (6.5 10-3 pc) X - X 14h,3 X BJJ\ 



(7) 



where V^h.a is the shock velocity in units of lO'^ km/s. 

Even though the actual size of the rims, A/Jj-imj may 
be larger than Ai?advi because of electron diffusion (see 
below), the constraint Ai?adv x P < Ai?obs is a simple 
consequence of the interpretation of the rims considered 
here. This straightforwardly translates into an upper limit 
on the synchrotron loss time and a lower limit on the 
magnetic field, _Bd, in the downstream region. Writing 
Ai?obs = ^obs,-2 X 10~^ pc, one obtains: 



T < 
' syn _ 



obs 



.5yr)x-xP-ixl/-i3xAi?obs,-2, (8) 



During the time Tgyn, the downstream electrons are ad- 
vected away from the shock at a velocity Vd = V^h/ ^7 '^syn is typically smaller than 3% of tlie age of tlie SNR, 



where P — P/4.6 is the reduced projection factor, and 
Sd > (210 fiG) X {APlrf/^ x V^^l x AR^^l^. (9) 
With the parameters in Tabled Eq. (|SJ| tells us that 
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^SNR, except for SN 1006, where the upper hmit is 6% (see 
Table 13) • This is of course consistent with the assumption 
that the maximum electron energy is not limited by ^snr, 
and it also justifies that we use throughout the current 
value of V^h- 

The lower limit on the downstream magnetic field will 
be made more precise below. Here, we simply use Eq. Q 
to derive a general limit related to the "advection thick- 
ness" of the rims, which we note Badv For Cas A, Kepler, 
Tycho, SN 1006 and G347.3 + 0.5, one obtains respec- 
tively (see also Table 01: 210, 170, 200, 57 and 61 /iG 
(for r = 4). This does not seem compatible with the 
mere shock compression of the upstream magnetic field, 
unless the latter is significantly higher than in the av- 
erage ISM and the compression ratio is of the order of 
7-10, as could resul t from non-linear CR shock modifica- 
tion l)Ellison et al. 1 12004'!. We thus argue that magnetic 
field amplification must occur around the shocks of young 
SNRs (see also below, and Volk et al., 2005, and references 
therein) . 

2.4. DifFusion downstream 

In addition to their advection with the downstream 
plasma, the energetic electrons diffuse in the local tur- 
bulent magnetic field, and if their effective diffusion ve- 
locity away from the shock is larger than Vd, they can 
fill regions larger than Ai?adv Writing again the sta- 
tionary version of the transport equation lO, but keep- 
ing the (uniform) diffusion term instead of the advection 
term, we obtain d^f/dx^ = f/DTsyn, and the solution 
/o(a;) cx exp(-x/Ai?diff), where Ai?diff = V -^Tsyn thus 
gives the diffusive scale of the rims. 

It was found above that the local magnetic field must 
have been amplified, therefore it is legitimate to assume 
SB » Bq, where Bq is the regular field, and we also as- 
sume a roughly isotropic turbulence. In this case, even 
though the value of the diffusion coefficient, D{E), is gen- 
erally not known, it cannot be smaller than the so-called 
Bohm value, Db — ri^v/3, as recalled in Scct.^ 

Db{E) ~ (3.3 X lO^^cm^s^i) x £;tcV x Bfo^f,. (10) 

With the help of Eq. one thus obtains: 

Ai?diff > (3.7 10-2 pc) X B-^T, (11) 

independently of the electron energy (as a result of the 
Bohm scaling). Just as in the previous section, we now 
obtain a lower limit on the downstream magnetic field by 
imposing the condition A_Rdiff x P < Ai?obs- This reads: 

Bd ^ (660 ^iG) X p2/3 X AR-^l'_,. (12) 

This result is interesting since it provides a lower limit 
on the value of i?d that can be derived directly from the 
observed thickness of the synchrotron rims, independently 
of the shock velocity. Moreover, since it was obtained 
with the Bohm diffusion coefficient, any other assumption 



SNR 
name 


max / I 

(x4P/r) 


Padv 

(mG) 


Pdifl 




Cas A 


< 2.6% 


210 


230 


1.1 


Kepler 


< 2.8% 


170 


180 


1.1 


Tycho 


< 2.1% 


200 


230 


1.2 


SN 1006 


< 5.9% 


57 


90 


2.0 


G347.3-0.5 


< 3.3% 


61 


77 


1.4 



Table 2. Constraints on the synchrotron loss time and the 
downstream magnetic field for the five young SNRs whose 
observed parameters are given in Table ^ In columns 3 
and 4, the constraints are lower limits derived from con- 
sideration of advection and diffusion only, respectively. To 
derive Pdiff we assume a Bohm diffusion coefficient, which 
leads to the less constraining result (to the lowest mag- 
netic field). In column 5, the superscript "B" indicates 
that the Bohm diffusion coefficient has also been assumed. 



would lead to a larger diffusion scale, and thus thicker 
X-ray rims. The Bohm diffusion coefficient therefore pro- 
vides the weakest possible constraint on the downstream 
magnetic field using the "diffusion thickness" of the rims, 
Pdiff, so that Eq. ((T^ is quite robust. Applying it to Cas 
A, Kepler, Tycho, SN 1006 and G347.3 + 0.5, with the 
values given in Table ^ one obtains, respectively: 230, 
180, 230, 90 and 77 ^J,G (cf. Table These lower hmits 
on Pd are remarkably close to those obtained with APadv 
and confirm that a strong amplification of the pre-existing 
field is required (Berezhko & Volk, 2004, and Volk et al., 
2005). 

2.5. Diffusion vs. advection 

While diffusion and advection have been considered sepa- 
rately above, it is interesting to note that the correspond- 
ing thicknesses of the rims produced by the highest en- 
ergy electrons accelerated at the shock should be of the 
same order of magnitude. Indeed, since APdig = Prgyn 
and APadv — VdTsyn, the condition APadv = APdig is 
equivalent to Tgyn = D/V^. Now D/V^ gives the order 
of magnitude of the acceleration timescale, Tacc, so that 
APadv and APdig are found to be comparable if and only 
if Tacc — ''"syn- This characterizes the electrons with the 
highest energy, Pc,max, in a loss-limited SNR (see also 
Vink, 2004). 

From Eqs. Q and (|ll|l . one easily obtains the thick- 
ness ratio in the case of a Bohm diffusion coefficient (at 
the observed gamma-ray energy of 5 keV, assumed to be 

close to Po.max): 

independent of the magnetic field. The numerical values 
obtained with our set of parameters are shown in Tabled 
for r = 4. Given the uncertainties in the measured quan- 
tities, they are remarkably close to l! This indicates that 
the value of the electron diffusion coefficient at Pe.max 



E. Parizot et al.: Observational constraints on energetic particle diffusion in young SNRs 



5 



cannot be significantly diS^erent from the assumed Bohm 
value. In the next section, we make this argument more 
precise by investigating the high-energy cutoff of the ac- 
celeration process, taking into account the time spent by 
the particles both upstream and downstream, and without 
restriction to the Bohm diffusion regime. 

3. Constraints on the SN acceleration parameters 

3.1. Acceleration timescale 

Apart from a few "lucky" particles that gain energy on 
a shorter timescale than the average in the stochastic ac- 
celeration process, the maximum energy, -£^c,maxi that an 
electron can reach is obtained by equating the synchrotron 
loss time to the acceleration time scale. According to the 
DSA theory, the latter is given by (Drury, 1983): 



3r frDd + Du 



- 1 V 



sh 



(14) 



where and are the upstream and downstream dif- 
fusion coefficients. 

Without any loss of generality, we can write = 
kuDs.u and I?d = kdDB,d, where D^^i is the Bohm dif- 
fusion coefficient in the local magnetic field (upstream or 
downstream), as given by Eq. (|10|l . and fcj is merely a 
number characterizing the difference between I?(-E'c,max) 
and the Bohm value: ki > 1. Note that in principle it de- 
pends on E, as in the case of Kolmogorov or Kraichnan 
type of turbulence, but this is not important here since we 
consider only the highest energy electrons. 

In order to limit the number of free parameters, we 
make the reasonable assumption that the upstream and 
downstream magnetic fields are related according to the 
usual jump conditions (shock compression of the compo- 
nents parallel to the shock front). Given the high mag- 
netic field values, which require field amplification, we as- 
sume that the magnetic turbulence is roughly isotropic 
upstream and write (following Berezhko, et al., 2002): 
Bd ~ v^(l -I- 2r2)/3Bu ~ 0.83rBu. Taking the down- 
stream magnetic field as a reference to compute the Bohm 
diffusion coefficient, one can rewrite Eq. p4|l as: 

(fcd + 0.83fcu) ^ 1.83/^0^^:^, (15) 



sn 



1 



vl 

sh 



where we have assumed, in the last equality, that ~ 
fed = ko, as expected if the amplified magnetic field has 
essentially the same structure upstream and downstream, 
apart from the above-mentioned shock compression. A 
more detailed investigation would involve additional free 
parameters describing the plasma properties on both sides 
of the shock and a complete model of turbulence, which 
is left for future work. Note, in particular, that a purely 
isotropic turbulence both upstream and downstream can- 
not hold in principle, because of the shock compression 
itself. It should also be noted that the effective value of 
the compression ratio "felt" by the particles depends in 



principle on their rigidity in the case of a strongly modi- 
fied shock. Only the protons of highest energy should feel 
the total compression ratio, rtot (including the precursor), 
while low-energy particles feel the much lower compression 
ratio of the sub-shock, r-gub, and the electrons of highest 
energy explore a region across the shock with an interme- 
diate effective compression ratio r^i (keeping in mind our 
basic assumption that the X-ray rim thickness is limited 
by synchrotron losses, so that these electrons have a lower 
energy than they would reach otherwise). We shall not go 
into such details here, as they would only complicate the 
formal treatment with no significant change in the results 
(e.g. always smaller than 30% for the deduced value of 
^'p.max! see below). In any case, we found that a situation 
with explicitly different values of the compression ratio for 
protons and electrons (and particles of different energies) 
in the case of a strongly modified shock results in an inter- 
mediate result between the test-particle case (no precursor 
nor sub-shock, and r = 4) and the pure r = 10 case. We 
shall thus confine our study below to such idealised cases, 
using only one ( "universal" ) compression ratio for all par- 
ticles. Numerically, Eq. 115|) gives: 



(30.6 yr) 



3r' 



16(r- 1) 



xko{E)xETovB^,lV^-^,. (16) 



3.2. X-ray cut-off and diffusion coefficient 

Likewise, to evaluate the effective synchrotron loss time, 
we need to distinguish the upstream and downstream 
energy losses (cf. Yamazaki et al., 2004). If Tu and 
are the time respectively spent upstream and down- 
stream, and Tsyn,u and Tgyn.d the corresponding syn- 
chrotron time scales, the effective energy decay law 
over one cycle, Tcycio — '''u + ^d, can be written 
^^(Tcycie) = -Eq exp(-ru/Tsy„^u) exp(-Td/rsy„,d) ■ The av- 
erage synchrotron loss time, (Tgyn), is obtained by identi- 
fication with i?(Tcyclo) = i?oexp(-rcyclo/ (7"syn)): 

(Tsyn) = (1.25 10^ yr) x iJ^eV X (i??oo)"' , (17) 

where (S^) = [B^t^ + BIta)/{tu + t^) is the average 
square field. With r,; cx Di/Vi, we have Tu/rd = D^/rDd ~ 
Bd/rB^, ~ 0.83 and 



{B') ^ B\ 



1 + l/0.83r2 
L83 



0.59 5^. 



(18) 



Thus, in terms of the downstream magnetic field: 



(T,y„) ~ (2.1 X 10^ yr) X Sfoo- (19) 
From Eqs. H16I) and H19|l . one obtains: 
i^cmax ^ (8.3 TeV) x f(r) x k'^'^ x B'^J^ F,h,3, (20) 



where f{r) = f{r)/f{A), with /(r) = v^r~~T/r: /(r) takes 
values between 1 and /(10)//(4) ~ 0.693. 

On the other hand, Eq. ^ gives -Ec,max as a function 
of the X-ray cutoff energy: 



E,, 



(21) 
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^0 (-E'Cjiiiax) 






r) in /iG 




-Be,max(Q,r) iu TeV 


a^in (Eq.OSl 




r = 4 


r = 10 


(1 ;4) 


(1 ; 10) 


(1/3 ; 4) 


(1/3 ; 10) 


(1 ;4) 


(1/3 ; 10) 


r = 10 


Gas A 


3.2 


1.5 


390 


280 


350 


250 


12 


15 


0.08 


Kepler 


4.5 


2.2 


340 


250 


300 


210 


11 


14 


0.09 


Tycho 


10 


4.9 


530 


400 


400 


300 


5.2 


6.9 


0.10 


SN 1006 


(0.40) 


(0.19) 


110 (84) 


95 (59) 


100 (82) 


91 (57) 


37 (32) 




0.37 


G347.3-0.5 


(0.87) 


(0.41) 


96 (93) 


84 (66) 


92 (89) 


79 ( 62) 


36 (37) 




0.13 



Table 3. Diffusion and acceleration parameters estimated for the five young SNRs considered in Tabled (see text). 



with the same dependence in B. Identifying both expres- 
sions, we can thus derive fco = -D(i?o,max)/-C'B(£'o,max) di- 
rectly from the SNR data, namely the X-ray cutoff energy 
and the shock velocity: 



= 0.14 X E. 



7,cut,kcV 



X Vl, X f{rr. (22) 



with 0.48 < /(r)2 < 1. The values obtained for r = 4 
and r = 10 are given in Table |21 kQ is found to be 
larger than 1 for Cas A, Kepler and Tycho, as required, 
or marginally lower than 1 for G347.3-0.5. In the case of 
SN 1006, fco ~ 0.4 for r = 4 and 0.2 for r = 10, which may 
favour low compression ratios and/or point to a possible 



overestimate of 



(e.g. if the synchrotron cutoff is 



less sharp than exponential) and/or underestimate of Vsh- 
Note that the cut-off energy and shock velocity given in 
Tabled were obtained from different regions in the SNR. 
The cut-off energies are quite uncertain, being mostly ob- 
tained by comparing X-ray and radio fluxes. They can 
be overestimated if the spectrum is concave (predicted 
by non- linear acceleration), or underestimated if the ra- 
dio flux is taken from a region larger than just the rim. 
Additional data on these parameters would be very valu- 
able for the present study. 

Despite the above uncertainties, it can be seen from 
Table |21 that fco keeps reasonably small values, which 
means that the actual diffusion coefficient at i?o,max is 
not very much larger than the Bohm value, and possibly 
very close to it in the case of SN 1006 and G347.3-0.5. Note 
also that numerical studies of diffusion in a turbulent mag- 
netic field found values of fco — 3~4 at the energy where 
rj, = Ar, the coherent length of a t urbulent magnetic field 
ijCasse et al. 1 120021: IParizot lEool . If the parameter val- 
ues for Cas A and Kepler are confirmed, this would favour 
a turbulent scale close to 7'L(£'o,max)- 



3.3. Rim thickness and diffusion regime 

We now repeat the calculations of Sect. |51 but taking 
into account advection and diffusion jointly, and using 
the above determination of the diffusion coefficient at 
Ec,ma^- To derive the scale of the high energy electron 
distribution downstream, we write again the stationary 
version of Eq. ((BJ, keeping advective and diffusive terms: 
Vdidf /dx) — Dd^f/dx^ — f /rgyn. The solution is again 
of the form f{x) cx exp(— ax), where a is the positive so- 
lution of the quadratic characteristic equation (Berezhko 



Volk, 2004): 



Da^ + Vda^ 




The scale of the emitting region is then Ai?rim = a ^ 
the observed (projected) size of the rims thus reads: 



obs 



P X 



2DIVa 



1 



and 



(24) 



By definition, this equation holds only at i?c,obs, given by 
Eq. 0, and we thus need to evaluate the synchrotron 
loss time and the diffusion coefficient at that energy. We 
already obtained Tsyn(-Be,obs) in Eq- ©• For the diffusion 
coefficient, we derived its value at i?c,max and now need 
make an assumption about the diffusion regime, i.e. the 
energy dependence of D. This is not known a priori, so we 
shall leave it unspecified and simply assume a power-law 
form: 



D{E)^D{Ea) ( — ) =koDB{E, 



E 



E, 



(25) 

where the index a is a free parameter. The Bohm regime 
corresponds to a = 1. Quasi-linear theory predicts a — 
1/3 (resp. 1/2) in the case of a Kolmogorov-like (resp. 
Kraichnan) spectrum of magnetic turbulence, and a sim- 
ple calculation leads to a = 2 for particles with a gyro- 
radius larger than the coherence length of the field. Both 
of these behaviors are very well reproduced with numer- 
ical simulations of particle diffusion (e.g. Parizot, 2004). 
Note that since D{E) must remain larger than Db{E), 
one must have fco(-E/i?o)"~^ > 1 at all energies where the 
power-law regime holds. Since fco ^ 10 at -Ec,max and since 
the injection energy in the acceleration process is many or- 
ders of magnitude below i?e,max, this implies that a < 1 
at least up to energies close to -Ec,max- 

Although neither i?c,obs nor -Ec.max are known yet, 
their ratio has a very simple expression. The typical syn- 
chrotron photon energy is simply proportional to the elec- 
tron energy squared, so by definition, with i5-y,obs — 5 keV: 



Ee^ohs 
Ecmi^.^ 



E. 



7, obs 



E. 



'7, cut 



1/2 



2.2 X E 



-1/2 
7,cut,kcV ■ 



(26) 



Since E^ cut < 5 keV for all five SNRs, we have E^ obs > 



E, 



(the ratios are respectively 2.0, 2.3, 4.1, 1.3 and 
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1.4). This can be used to set a lower limit to a, by requir- 
ing that D{Ec^ohs) > £'Bohm(-Bc,obs)- From Eq. this 
translates into: 



a > 1 - 



lnfco(£^c^inax) 
ln(i?o,obs/£'c,max) 



(27) 



This is however not constraining in the case of Cas A, 
Kepler and Tycho, since the lower limit is negative. For SN 
1006 and G347.3-0.5, we already noted that the value of kg 
derived from the data is lower than 1. In the following, we 
shall assume ko = 1 in this case, putting the unrealistic 
computed value down to reasonable uncertainties in the 
measured parameters. 

3.4. Self-consistent magnetic field and maximum 
electron energy 

Replacing Eq.^^ in and using we get: 



D{E,^obs) = (1.0 10^4 cm2s-i)S~'o/V,2_3i^(a,r), (28) 

where we defined 

Kia, r) = f{rf x 2.2" x E-^J^,^^^ (29) 

Using Eq. I^), the ratio in the denominator of Eq. H24|l 
follows: 



4J(^e,obs) 

,obs ) 



8.l(0V(a,r), 



(30) 



which is independent of both and V^h- Replacing 
Eqs. H30I) and H28I) in H24I) . we obtain a self-consistent ex- 
pression of the downstream magnetic field as a function 
of a, r and the measured SNR parameters: 



Bd ~ (520 /LtG) 



{r/A)PK{a,r)V,^^:i/\R~l_^ 
v/l + 8.1(r/4)2ii:(a,r)- 1 



2/3 



(31) 



The results are given in Table 13 for different compression 
ratios and diffusion regimes. For SN 1006 and G347.3-0.5, 
we assumed fco(-E'o.max) = 1 (i-e. D = I?Bohm) and gave the 
value obtained with the (non physical) computed value of 
fco < 1 between parentheses. As already noted, the esti- 
mated magnetic fields are much larger than those typi- 
cal of the interstellar medium and require strong ampli- 
fication at the shock. Lower magnetic fields are obtained 
with lower values of a and/or larger compression ratios. 
Two limiting cases are considered; (a — \,r ~ 4) and 
(a = 1/3, r = 10). 

With the above value of the magnetic field, we can now 
derive the maximum electron energy in a fully consistent 
way, from Eq. if^ : 



E, 



(9.6 TeV) X E, 



1/2 

7,cut,kcV 



(r/4)Pg(a,r)Kh.3Aflobs.-2 
v/1 + 8.1(r/4)2i^(a,r) - 1 



-1/3 



(32) 



The results are given in Table El they are of the order 
of 10 TeV for Cas A and Kepler, 5 TeV for Tycho, and 



40 TeV for SN 1006 and G347.3-0.5, almost insensitive to 
the parameters a and r (the two extreme cases are shown). 

Finally, the value of £^e,max and the associated value 
of -Ee.obs allow us to set a lower limit to a, by requir- 
ing that the total acceleration time of electrons up to the 
highest energies do not exceed the age of the SNR, tsNR- 
Integrating the acceleration timescale, Eq. H14() . from Ei^j 
to i^cobs > ^^inj, using Eq. we find (for a > 0) 



^acc {-f^i 



c,obsJ 



Bin 



dE 1 

Tacc{E)—— ~ — Tacc(^^c,obs)- (33) 

E a 



Now from Eqs. and we get: 

racc(^^e,obs) ^ {93 jr) f'^ B^^^^^ K {a, r) 

c.(93yr)Br„f x2.2"i^;J/^r^^ 

so that the condition iacc(-£'c,obs) < ^snr reads: 



(34) 



apP{r/A) ^ 7.8 yr AEpbs.-z ^3^.^ 



y/l + 8.1{r/A)^K{a,r) - 1 tsNR Vs 



sh,3 



Solving for a, we get the lower limit amin- For r = 4, this 
is of the order of 0.05, which is not really constraining. 
The values of ckmin are given in Tabled for r ~ 10. They 
also are not constraining, and compatible with either a 
Kolmogorov-like or a Bohm diffusion regime (note that 
the case r = 10 with a < 1 is not relevant for SN 1006 
and G347.3-0.5 anyway). 

4. Maximum energy of the accelerated protons 

In the previous section, we estimated the value of the am- 
plified magnetic field in the vicinity of the shock and of the 
diffusion coefficient at the highest electron energy, -Ec.max- 
We now consider their implications for the highest proton 
energy, -Ep^max, which (in the most favourable case) is lim- 
ited by the age of the SNR and is thus obtained from the 
condition iacc(-£'p,max) = ^SNR- Integrating like in Eq. H33|l 
up to i?p,max, one obtains: 



^acc (-Epjinax) 



'^acc (Eo) (E, 



and thus 



E,, 



Eo 



a^SNR 



^p,max 
En 



1/a 



(36) 



(37) 



ciEo). 

Choosing Eq — E'cmax or i?o,obs and replacing the cor- 
responding expressions, we derive ii^p,max as a function 
of the only remaining free parameters, a and r (through 

-Bioo): 



Er, 



(49TeV)Sroy' 



tSNR OlB 



3 /2 

100 ^l/2+a/2 



93 yr 2.2^ 



7,cut,kcV 



1/a 



(38) 

where B is given by Eq. (|31|l . This maximum proton en- 
ergy is shown as a function of a for the SNRs under 
consideration in Fig. The shape of the curve is eas- 
ily understood by considering the interplay between the 
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alpha alpha 

Fig. 1. Maximum proton energy, i?p,maxj as a function of the diffusion regime index a, for r = 4 (left) and r = 10 
(right). The dashed part of the curves correspond to values of a that lead to a diffusion coefficient smaher than the 
Bohm limit at -Ep^max, and are therefore excluded. The highest possible value of iSp.max is materialised by the circles 
giving the value ckopt above which diffusion regime indices are valid - notably the Bohm index a = 1, but lead to 
smaller values of E'p^max- 



SNRs 


r = 4 


r = 10 


modified SNR parameters 


Otopt 
{— Omin) 


-Bp, max (in PeV) 


Qopt 

(= Qmin) 


Sp,niax (in PeV) 


Qopt 

(= Omin) 


Sp,max (in PeV) 


a = "opt 


a = l 


a = aopt 


a = l 


a = Oopt 


a = l 


Gas A 


0.72 


0.75 


0.35 


0.87 


0.32 


0.25 


0.71 


2.1 


0.64 


Kepler 


0.64 


0.77 


0.29 


0.76 


0.33 


0.21 


0.65 


2.2 


0.53 


Tycho 


0.52 


0.63 


0.15 


0.58 


0.27 


0.11 


0.56 


2.0 


0.28 


SN 1006 


1 


0.29 


0.29 


1 


0.13 


0.13 


1 


0.73 


0.73 


G347.3-0.5 


1 


0.82 


0.82 


1 


0.34 


0.34 


0.94 


2.0 


1.7 



Table 4. Derived values of a and i?p.max for the five SNRs considered in Tabled (see text). In addition to the cases 
of a compression ratio of 4 and 10, we consider an "extreme case" with SNR parameters artificially modified so as 
to increase the value of £'p,max: this amounts to increasing the shock velocity by 25% and dividing the observed rim 
thickness by a factor of 2. In each case, we give of i?p,max for aopt and for a = 1 (i.e. Bohm diffusion regime). 



instantaneous acceleration rate, Tacc, and the integrated 
acceleration time, iacc, at -Ep,max- Above, we were able to 
calculate the diffusion coefficient at i?o,inax, which gives a 
fulcrum for D at all other energies. For very low values of 
a, the diffusion coefficient at low energy is much higher 
than with the Bohm diffusion, and the integrated accel- 
eration time becomes very long. As a increases, this situ- 
ation gets better, but the ratio between r'(ii^p,max) and 
£'(£^o,max) also incrcascs, so that Tacc(£^p_niax) becomes 
higher and it takes longer to reach higher energies. The 
best compromise is for a ~ 0.1-0.2. However, these val- 
ues of a are excluded because they lead to diffusion co- 
efficients at -Bp, max that are much lower than the Bohm 
diffusion coefficient. Therefore, the curves in Fig. ^ are 
only valid for pairs of values of a and -Ep ^ax such that 
£'(-Bp,max) > £'Bohm(£^p,max)- Exactly as in Sect. ESI we 
must make sure that fco(-Ep_i„ax/-E'o,max)"~"'^ > 1, which 
corresponds to the part of the curves shown in plain 



line (dashed line otherwise). Therefore, although rela- 
tively high values of -Ep^max would seem to be reached for 
a ~ 0.1-0.2, notably higher than the energy of the knee, 
these values are not physical and must be excluded. 

For each SNR, one may define the "optimal value" of 
a as the value aopt giving the highest possible value of 
-Ep.max- These correspond to the smallest "allowed" val- 
ues of a as discussed above, and are given in Table 0] 
together with the highest proton energy obtained under 
the assumption of a Bohm regime. As can be seen, aopt is 
typically between 0.5 and 0.7 for r = 4 and between 0.6 
and 0.9 for r = 10, except for SN 1006 and G347.3-0.5, 
where we already noted that a < 1 is not favoured. As 
can be seen, even in the optimal case, the highest possible 
proton energy falls short of the energy of the knee in the 
cosmic-ray spectrum, at 3 10^^ eV, by about a factor of 4 
(or even 10 in the case of SN 1006). This is even more true 
if efficient acceleration is assumed, so that non-linear ac- 
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100 rj, 



10 



> 



a. 




Fig. 2. Same as Fig. ^3-, with artificially modified SNR 
parameters: V^h higher by 25% and Ai?i.ii„ two times lower. 



celeration effects lead to a higher compression. For r = 10, 
-E-p^niax is systematically lower than the knee energy by one 
order of magnitude. It should also be noted that a Bohm 
diffusion regime would give lower values of -Ep.max, which 
is due to the fact that -©(E'p^max) is then larger than the 
Bohm value. The diffusion regime index that gives the 
lowest possible diffusion coefficient at i?p,max is precisely 
Qfopti hence the higher maximum energy. 

From the above analyses alone, our results do not set 
any interesting constraint on the actual value of a. Values 
lower than aopt may be considered as perfectly valid and 
merely resulting in lower values of the maximum proton 
energy. For instance, one might imagine a situation where 
electrons and protons have the same maximum energy, 
with the maximum turbulence scale fixed by i5c,max so 
that protons with a higher energy are not scattered res- 
onantly and are thus not accelerated efficiently. However, 
we now claim that a further step can be taken if we con- 
sider the underlying turbulence generation process. While 
we just mentioned the possibility that the maximum scale 
of turbulence be set by electrons alone, so that the pro- 
tons cannot be accelerated above i?c,max, there is no reason 
a priori why such a situation would occur. If the ultra- 
relativistic electrons are able to generate the turbulence 
necessary to keep the acceleration process going, then so 
should the protons, since at these energies they are essen- 
tially equivalent from the electromagnetic point of view. 
Besides, electrons only carry a negligible fraction of the 
energy of relativistic particles. It is thus reasonable to 
assume that whatever underlying processes actually lead 
protons to £'c,max, they should not cease to work at that 
particular energy, which is determined by limiting pro- 
cesses affecting electrons only, nor at any other energy 
bearing no physical meaning for the process under con- 
sideration. In other words, the only natural limitation for 
the proton maximum energy is related to the acceleration 



time, tacc(^'p) < ^SNR (since the size-related limitation is 
less stringent for the SNRs considered here). 

Applying this argument to the present study, we can 
actually reject any value of a lower than aopt, since for 
these values the maximum proton energy is not limited 
by the age of the SNR, but "artificially" by imposing a 
cutoff in the turbulence spectrum at the energy where 
D{E) = r'Bohm(-E'). If that energy has to have a physi- 
cal meaning, in terms of the Larmor radius of the highest 
energy protons (limited by the SNR age) in the amplified 
magnetic field, then a must be larger than the very value 
that we noted aopt above. From this point of view, i5p,max 
is known to be given by Eq. H37() . and the requirement 
that L'(£^p,max) be larger than £'Bohm(£^p,inax), together 
with the observational constraint on £'(i?e,max), straight- 
forwardly leads to a > aopt, i-e. to an interpretation of 
aopt as a lower limit for the diffusion coefficient index (as- 
suming that the index does not significantly depend on 
energy, which may be expected to be the case up to en- 
ergies close to the maximum resonant energy, such that 
tl — Ainax; Casse et al., 2002). Higher values of a remain 
possible in principle, since there is no reason a priori to 
require that the diffusion coefficient at the highest proton 
energy (limited by tsNR.) have exactly the Bohm value. 

Another interesting result is that i?p,max is roughly 
identical for all the SNRs considered, with values around 
0.7 PeV at most, despite the different values of their pa- 
rameters. This may indicate that not only these SNRs, 
but most of them probably cannot accelerate protons up 
to the knee. Since all the processes involved only depend 
on the rigidity of the nuclei, heavier nuclei of charge Z 
should reach energies Z times larger, i.e. ^ 20 PeV for Fe 
nuclei in the optimal case. 

In order to investigate the errors involved in our calcu- 
lation due to uncertainties on the measured parameters of 
the SNRs, we considered what may be thought of as an ex- 
treme case, where the parameters have been moved consis- 
tently towards values leading to higher maximum proton 
energies. In Fig.|21 we show i?p,max(a) for this case where 
we have increased the shock velocity by 25% and divided 
by 2 the rim thickness, keeping the most favourable value 
of the compression ratio, namely r = 4. The influence of 
the X-ray cut-off energy is somewhat more complicated, 
since a higher value of Ey^cut results in an upper shift of 
the curves in Fig. ^ thereby increasing -Ep.max at a given 
a, but at the same time the allowable range of values of a 
is reduced, and aopt is shifted upwards, reducing the max- 
imum E'p.max- The net effects of an upper (resp. lower) 
shift of i?-y,cut by a factor of 2 is then an increase (resp. 
decrease) of i?p,max by less than 10%. Since the adopted 
values of Ej,cut may already be considered upper limits 
(see above), we did not change this parameter in Fig. 13 
The corresponding values of i?p,max for "opt and for the 
Bohm diffusion regime are given in Table ^ As can be 
seen, even in this case i:^p,max(ctopt) remains consistently 
below the knee, although only by 40% for all SNRs except 
SN 1006. For a = 1, £;p,max is stiff a factor of 4-10 below 
the knee for all SNRs except G347.3-0.5, where the factor 
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is ~ 2. In other words, it appears that the strong magnetic 
field amplification reported here is still not enough to en- 
able young, isolated SNRs as considered here to provide 
(at least directly) the main contribution to the observed 
Galactic cosmic-rays above 10^^ eV, i.e. in the knee region 
and beyond. 

5. Conclusion 

In this paper, we have investigated which constraints on 
the parameters of young SNRs could be derived from the 
observed thickness of their X-ray rims, under the assump- 
tion that this thickness is limited by the synchrotron losses 
of the highest energy electrons - which is fully consistent 
with our results. This extends the analysis of Berezhko & 
Volk (2004) simply by relaxing the assumption of Bohm 
diffusion, and then provides an opportunity to constrain 
the diffusion coefficients from the data, in addition to 
other SNR parameters. Considering both the diffusion and 
advection of electrons in the downstream medium, we have 
shown that the magnetic field must be amplified in all the 
SNRs considered, with downstream values between 250 
and 500 /uG in the case of Gas A, Kepler, and Tycho, as- 
suming a shock compression ratio of 4 (unmodified strong 
shock). If one assumes a compression ratio as high as 10, as 
may be the case if substantial shock modification occurs 
(very efficient shock acceleration), magnetic field values 
about 25% lower arc obtained, still requiring strong field 
amplification at the shock. 

The value of the diffusion coefficient at the highest 
electron energy, ii'e.max, can also be derived from the data, 
by relating the X-ray energy cutoff to the acceleration 
timescale. We found that D{Ec,max) is typically between 
1 and 10 times the Bohm diffusion coefRcicnt. We were also 
able to derive the value of i?c,inax, which ranges between 
5 and 50 TeV, roughly independently of the compression 
ratio. 

This analysis essentially constrains the diffusion co- 
efficient at -Ec,max and provides no significant informa- 
tion about the diffusion regime, i.e. the dependence of 
D on energy. The comparison of the electron acceler- 
ation timescale with the age of the SNR showed that 
the data is compatible with any value of the energy de- 
pendence power-law index in the most interesting range: 
1/3 < a < 1. However, by considering in the same frame- 
work the acceleration of protons and notably the maxi- 
mum proton energy, we were able to define an "optimal 
value" of a, typically in the range 0.5-1., for which the 
highest proton maximum energy could be reached. Further 
noting that the only natural limitation to £'p,max is the age 
of the SNR (in the cases under consideration), we also ar- 
gued that the above value of aopt actually sets a lower 
limit to the possible diffusion coefficient index in these 
SNRs. In other words, values of a lower than aopt would 
lead to diffusion coefficients lower than the Bohm value 
at E'p.max, as determined by tacc(-Ep,max) = isNR- This 
favours diffusion regimes between the Kraichnan regime 



{a = 1/2) and the Bohm regime (a = 1), rather than a 
Kolmogorov-like regime, with a = 1/3. 

Finally, we deduced the maximum energy of nuclei ac- 
celerated in the five young SNRs under study. Although 
this energy depends on the diffusion regime, a general up- 
per limit could be derived. We found E^^x — Z x 0.7 - 
0.8 PeV for all SNRs (except for SN 1006: 0.3 PeV), ob- 
tained with the optimal value of a 0.5-0.7 (or 1 for SN 
1006 and G347.3-0.5). Assuming a Bohm diffusion regime 
would lead to lower values of -Ep,max, around 0.3 PeV. 
Likewise, a higher compression ratio (r = 10) would cor- 
respond to ii'max ^ Z X 0.3 PeV in the optimal case. 

Of course, the values calculated here depend on the 
accuracy of the observed SNR parameters, namely the 
synchrotron cut-off energy, the velocity of the expanding 
shock and the size of the X-ray rims. We investigated the 
influence of a misestimate of these parameters (the most 
influential of which is -B-y.cut) by quite large amounts and 
found only small modiflcations of the highest proton en- 
ergy, remaining marginally, but consistently lower than 
the knee energy. It should be remembered, however, that 
all our results were obtained under the assumption of a 
uniform, isotropic turbulence. We nevertheless point out 
that non uniform turbulence with the thickness of the rims 
limited by the relaxation scale of the turbulence (Pohl et 
al. 2005) would actually relax the constraint on the down- 
stream magnetic field and make the acceleration less effi- 
cient, especially at large gyroradii, leading to even lower 

values of -Bp^max- 

In conclusion, the derived limit on the maximum en- 
ergy of protons and nuclei accelerated in the young SNRs 
considered here appears to be quite robust, of the order 
of 800 TcV or up to 2 PeV by combining all uncertainties 
towards that direction, even though the magnetic field has 
been considerably amplified at the shock. In the most pop- 
ular framework where the Galactic component of cosmic- 
rays extends up to the ankle, i.e. up to ~ 3000 PeV, our 
results would thus suggest that if isolated SNRs are the 
main sources of GCRs, an additional GR component is 
required above the knee(s). More specifically, it appears 
very difficult for the considered SNRs in their current 
stage of evolution to produce protons up to the knee of 
the cosmic-ray spectrum, at ~ 3 10^^ eV, and essentially 
impossible to accelerate Fe nuclei up to either the ankle 
or the second knee at ^ 500 PcV. However, if the knee 
is associated with the contribution of a single nearby su- 
pernova (Erlykin & Wolfendale 1997, 2004) dominating 
the bulk of the Galactic cosmic-rays produced by some 
other mechanism, the expected contribution, as derived 
here directly from measured parameters, would show the 
successive cut-off of the most abundant nuclei according 
to E^axii) c^iZiX 0.5-1 PeV. 
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